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We demonstrate that homogeneous naturally-occurring materials can form non-magnetic negative refractive 
index systems, and present specific realizations of the proposed approach for the THz and far-IR, frequencies. 
The proposed structure operates away from resonance, thereby promising the capacity for low-loss devices. © 
2012 Optical Society of America 



Following the initial proposal by Veselago in 1968J2 
negative refraction materials spent over 30 years as a 
forlorn curiosity before being resurrected with renewed 
interest from both theoretical and experimental groups. 
Within the last several years it was realized that these 
materials (known also as left-handed materials) possess 
unusual properties, some of which were not recognized at 
the time of their conceptions. These properties include 
resonant enhancement of evanescent fields, potentially 
enabling near-perfect imaging below the diffraction limit 
and leading to a new class of optical devices,^ as well as 
nontrivial behavior in the nonlinear regimeP Despite ini- 
tial controversy over the realizability of negative index 
materials (NIMs), successful proof of principle demon- 
strations have been accomplished HMHHH3 

Existing designs for left-handed materials rely 
on achieving overlapping dipolar and magnetic res- 
onancespOGO or using photonic crystals near the 
bandgapJ-!30 Both of these approaches necessitate com- 
plicated 3D patterning of the medium with microstruc- 
tured periodic arrays. Fabrication of such structures 
presents significant challenges even for GHz applications, 
while manufacturing systems for higher frequencies be- 
comes harder still. Furthermore, near-resonant opera- 
tional losses impose severe limitations on the imaging 
resolution.^ 

As an alternative to periodic systems, a waveguide- 
based implementation of a NIM was proposedpl which 
obviates the need for negative magnetic permeability and 
does not require periodic patterning. This approach cir- 
cumvents major manufacturing obstacles to achieving 
NIM behavior at terahertz or optical frequencies. 

To achieve this behavior, the waveguide material must 
possess characteristics of a uniaxial medium with a sig- 
nificant anisotropy. Furthermore, this anisotropy must 
ensure that e± (the component of e transverse to the 
planar waveguide) is negative, while eii (in-plane com- 
ponent) remains positive. TM modes in such waveguide 
undergo negative refraction in the waveguide plane, and 
propagate with negative phase velocity.^ 

One of the key aspects in designing such an NIM 
system is selecting the material for the waveguide's 
anisotropic core. Several options have been proposed for 



the core material, in particular, nanostructured compos- 
ites in a dielectric host and quantum well structures.^ 
While being within the grasp of existing technology, the 
fabrication of such systems remains highly challenging. 

In this Letter we present an alternative approach to 
non-magnetic NIMs for THz and far-infrared domains 
based on naturally occurring materials with large di- 
electric anisotropy. In particular, we discuss the possi- 
bility of negative refraction in a waveguide-based system 
at approximately 20 /im, 58 /im and 255 fim using, re- 
spectively, sapphire, bismuth, or triglycine sulfate in the 
waveguide core. We focus on monocrystalline bismuth as 
an attractive option for manufacturing the NIM waveg- 
uide core thanks to its large anisotropy and availability 
of samples with high purity. 

In a planar waveguide with anisotropic dielectric core 
the wave vector components k z and k y are governed by 
the dispersion relation 
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with v = (l — /t 2 c 2 /eii a; 2 ) where uj is the frequency of 
light, e = ex (en) for the TM(TE) modes, k is the trans- 
verse mode parameter, and k z and k y lie in the waveguide 
plane. For perfectly conducting waveguide walls (a good 
approximation for silver and other metals at THz and 
far-IR frequencies^), k = mir/d, where m is an integer 
and d is the thickness of the waveguide.^ The effective 
refractive index for propagating waveguide modes in this 
system is given by n 2 s = e v. To support propagating 
modes, e and v must have the same sign. The case e > 0, 
v > is typically realized in an isotropic planar waveg- 
uide operating above cut-off. However, in the case e < 0, 
v < negative refraction occurs,^ with refractive index 
given by 

n cf f = -y/ex v. (2) 

Note that if k is regarded as the transverse wave vector 
component, Eq.pJ) can be rewritten as a hyperbolic dis- 
persion relation fc^/e 1 1 — k 2 / \ e±\ = lu 2 /c 2 , identified in 
Ref. ~ 



13 with the onset of negative refraction. 
As follows from the definition of v, 



the propagation 
with negative phase velocity occurs only for TM modes 
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Fig. 1. (A): The real (top panel) and the imaginary (bottom panel) parts of the dielectric function of TGS; the 
monoclininc C2 axis is along the "perpendicular" (_L) direction. (B): Same for sapphire; the crystallographic c axis 
is along the "perpendicular" (_L) direction. (C): Same for monocrystalline bismuth. 



when eii > 0, ej_ < 0. This behavior is observed in a num- 
ber of materials where structural anisotropy strongly af- 
fects the dielectric response. 

One example of such materials is triglycine sulfate 
(TGS), a compound widely used in fabricating infrared 
photodetectors. Spectroscopic studies of the crystal at 
low temperature have shown that phonon modes polar- 
ized parallel to the crystal's monoclinic Ci axis signifi- 
cantly differ in frequency from phonons transverse to the 
axis. This results in a large anisotropy in the dielectric 
tensor along these directions. In particular, dielectric re- 
sponse for the field polarized along the C2 axis features 
a resonance at 268 /xm, which is absent if the incident 
field is polarized transverse to the C2 axisP^ Dielectric 
function ej_ in the vicinity of this resonance can be fit- 
ted with the Lorentz-Drude modelp^ while e± in this re- 
gion can be taken approximately constant EMU] Lorentz- 



Drude model parameters from Ref. 15 were used to con- 
struct Fig. [lja). As is evident from the figure, e±_ < 0, 
while en > in the region 250 < A < 268 /im. Fur- 
thermore, the imaginary part of e becomes small away 
from the resonance, minimizing absorption. A TGS-filled 
waveguide with C2 axis oriented perpendicular to the 
waveguide plane would support negative index propa- 
gation, while suffering from minimal propagation losses 
(Im[e]~l(T 3 at 250 /mi). 

Whereas the phonon anisotropy of TGS exists in 
the low-THz domain, for other materials, it may occur 
in a different spectral band. In particular, the strong 
anisotropy of the dielectric response of sapphire (AI2O3) 
is also due to excitation of different phonon modes (po- 
larized either parallel or perpendicular to the c axis of 
the rhombohedral structure), but occurs around 20 /xm. 
Fig. [l|b) shows experimentally-determinecP^ e\\ and e± 
as functions of frequency. As with TGS, a region of e± < 
0, €|| > is evident in the experimental data. This poten- 
tially enables a sapphire-based waveguide NIM (with the 
c axis of sapphire core perpendicular to the waveguide 
plane). Note that the minimum of the material absorp- 
tion occurs in the frequency range of interest. 



Anisotropic phonon excitations are not the only mech- 
anism that can lead to strong dielectric anisotropy. Bis- 
muth, a Group V semimetal with rhombohedral lattice 
and trigonal symmetry, exhibits such anisotropy due to 
a substantial difference in its electron effective masses 
along different directions in the crystal. 

In the frequency region of interest, the spectral de- 
pendence of the electric permittivity of Bismuth can be 
adequately described by the Drude model, 



e L 1 



UJ Z + lUJT 
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with the lattice permittivity, oj p \ — Ne 2 /e^iries the 
plasma frequency, and r the relaxation time. These pa- 
rameters are known from interferometric and reflectance 
studies of Bi samples. In particular, plasma frequency of 
pure Bismuth at 4 K was measured to be 158 cm -1 for 
the incident iS-field polarized perpendicular to the trigo- 
nal axis, and 186 cm' 1 for the field polarized parallel to 
the axisP^ These values are in agreement with other ex- 
perimentsliHII22l The lattice dielectric constant for the 
field perpendicular to the trigonal axis was found to be 
110±10 cm -1 j22l in reasonable agreement with Ref. \l9 



For polarization parallel to the trigonal axis, £^=76.^ 
There can be substantial variation in the relaxation 
time t depending on the purity of the sample. We take 
t=0.1 nspSl however, this is a conservative estimate; for 
low temperatures, relaxation times over an order of mag- 
nitude greater have been reported as far back as 1975.^ 
Even with t=0.1 ns, the typical ratio of imaginary and 
real parts of the dielectric function in Bi is on the or- 
der of 0.1% in the frequency interval of interest, which 
enables many imaging and transmission applications.^! 
It should also be noted that high-quality single- crystal 
films as thin as 1 /im, with the trigonal axis (C3) oriented 
perpendicular to the film plane, have been reported,^ 
thereby essentially solving the technological issues in fab- 
ricating the proposed negative index device. 

Fig. [ljc) shows the behavior of real and imaginary 
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Fig. 2. The effective refractive index of the three lowest 
order modes in a bismuth waveguide. Inset: numerical 
simulation showing the refraction of a beam within a 
metallic waveguide (of thickness d — 4.5 (im) at an in- 
terface between an isotropic dielectric with e = 55 and 
monocrystalline bismuth; A — 61 /im. 



components of e for Bi based on Eq.([3]). The most promi- 
nent feature of these plots, the transition from e > to 
e < 0, is determined by the highly anisotropic plasma 
frequency. This anisotropy creates a window between 
A = 53.7 /im and 63.2 /im where e < for the -E-field 
along the C3 axis, while e > for E transverse to C3. 
The existence of such 10 /im window was confirmed by 
direct measurement P3 

To allow for left-handed propagation in this frequency 
interval, Bi should be integrated into the core of a pla- 
nar waveguide, with the C3 axis oriented in the trans- 
verse direction. In Fig. [2] we examine the behavior of the 
effective refractive index n c g [Eq.Q] for the proposed 
sub-critical (d < A/2) waveguide structure with bismuth 
core. Note that negative effective index is possible for all 
modes over the entire (e_i_ < , en > 0) range. Negative 
refraction behavior of our system was further confirmed 
by a numerical calculation of the electric field incident 
on the Bi waveguide. The results of this calculation are 
presented in Fig. [2] (inset). We assume a TM wave with 
a gaussian profile, mode-matched into the Bi waveguide 
in the transverse direction by e.g. propagating the beam 
from a metallic waveguide of the same thickness, filled 
with a regular dielectric. One can clearly see the nega- 
tive refraction at the boundary. Furthermore, it is evi- 
dent that attenuation of the transmitted wave is weak, 
as expected from low values of the imaginary part of the 
dielectric constant [Fig.[ljc)]. 

In addition, this calculation shows that the NIM 
waveguide remains transparent despite the fact that 
transverse dimension of the waveguide is much smaller 
than the wavelength, which indicates strong confinement 
of the field within the core (since the cladding is assumed 
to be perfectly conducting). This behavior is not found 
in a subwavelength dielectric waveguide, where much of 
the field spreads into the cladding, or a subwavelength 
metallic waveguide, which does not support propagating 
modes. Such strong field confinement may find applica- 



tions in photonic structures and nonlinear opticsP^l 

In conclusion, we have proposed a novel negative re- 
fraction system for several wavelengths from low-THz to 
far-IR. Our approach is non-magnetic, avoids the use of 
periodic patterning, utilizes naturally occurring materi- 
als, and promises the capacity for low-loss devices. 
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